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TECHNICAL NOTE 3686 


EXPERIMENTAL MEASUREMENTS OF FORCES AND MOMENTS ON A 
TWO-DIMENSIONAL OSCILLATING WING 
AT SUBSONIC SPEEDS 1 

By Sherman A. Clevenson and Edward Widmayer, Jr. 


SUMMARY 


Experimental results of lifts and moments about the quarter chord 
of a two-dimensional wing at subsonic Mach numbers are presented. A 
comparison of the experimental magnitude of the lift vector with the 
theory as given by Dietze showed good agreement. Comparisons of the 
moments and the quadrature component of lift with theory indicated that 
some refinements in the testing technique are necessary for the experi- 
mental determination of these quantities in the transonic range. 


INTRODUCTION 


An extensive amount of work has been done on the theory of air 
forces on a harmonically vibrating wing in a two-dimensional compres- 
sible medium by Possio, Frazer and Skan, Dietze (references 1 to 4), and 
others. Flutter of wings in a compressible medium has also been treated 
by Garrick (reference 5)j Frazer and. Skan (reference 2 ), and others in 
which the aforementioned theoretical air forces are used. Numerous 
wind-tunnel tests have also been made on the flutter of cantilever wings 
in the subsonic speed range (reference 6 ). These tests show reasonable 
agreement with the indicated theoretical trends. Some work on the 
experimental determination of the oscillating air forces at high 
subsonic speeds has been done by Bratt and Chinneck (reference 7 ), and 
the oscillating moments on a small airfoil oscillating about the mid- 
chord position were determined at high subsonic speeds for a range of 
reduced frequencies of 0.0025 to 0.015. The aerody nami c coefficients 
have been calculated by Dietze for a reduced-frequency range to cDj. = 0.6 

and for a Mach number range to M = 0 . 7 . 

The present paper deals wit'h some experiments at subsonic speeds 
in which a two-dimensional wing was oscillated about the quarter- 


-^Supersedes declassified NACA Research Memorandum L9K28a by 
S. A. Clevenson and E. Widmayer, Jr., 1950* 
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chord position in a reduced frequency range from 0.15 to 0.35* Since 
extension of the theory into the transonic range is extremely difficult 
and open to considerable douht because of the occurrence of shocks and 
mixed flow on the wing, it appears desirable to perform an exper ime nt 
for determ in ing the oscillating air forces and moments in this region. 
Although this paper presents few data in the transonic region, it does 
present an experimental method for obtaining these forces and moments. 

A comparison is made of the experimental forces and moments with the 
theory of Dietze. Although comparisons could not be made for all 
components, it is believed that the presented method can be refined to 
give sufficient accuracy for the determination of the components of lift 
and moment of an oscillating airfoil in the transonic range. 
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semichord 
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mass density of gas 

velocity of gas 

unit span of wing 

oscillating lift (L 0 sin(asb + q>p)) 

oscillating moment about quarter chord (Q 0 sin(ajfc- + cpg)) 

Mach number 

time 

amplitude of pitching oscillation (oq sin ajfc) 
natural circular bending frequency 

natural circular torsion frequency 

circular frequency of oscillation 
reduced frequency (bco/v) 

slope of lift curve 
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The wing has a span of 2b Inches and a chord of 8 inches, with an 
NACA 63— 010 airfoil section. It was of fabricated construction, having 
a steel spar with duralumin skin and was mass balanced about the 
quarter chord to eliminate any mass— inertia force. This airfoil was 
supported on bearings at each end and had only freedom to pitch about 
the quarter-chord point. The wing oscillated about a mean angle of 
attack of 0° and had a maximum amplitude of ±2.30°. Two variable- 
frequency, electromagnetic shakers phased to produce pure moment an d 
a torsion rod were attached near one end of the wing (fig. 1). The 
other end of the torsion rod was fastened to the tunnel wall. Various 
reso n an t frequencies could be obtained by the use of torsion rods 
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haring different torsional stiffness. The frequencies of the wing— rod 
combinations used are as follows: 


Wing— rod combination 

“h 

(radians/sec) 

<jd„ 

a 

(radians/sec) 

1 

528 

194 

2 (Data not reported) 

546 

275 

3 (Data not reported) 

572 

361 


Tunnel 

The test section of the Langley 4^-— foot flutter research t unn el 

was modified hy means of auxiliary walls so that the test section was 

24 inches wide. A — r — inch gap existed between the model and the tunnel 

16 

walls. All tests reported were made in Fre ore-12 in which the speed of 
sound was about 520 feet per second. The use of this test medium 
permits the attainment of more than twice the reduced frequency as 
obtained in air for a given Mach number and torsional frequency. 

Reynolds numbers for these tests were from 1 x 10 ^ to 5.5 x 10^. 


Instrumental; ion 

The instrumentation may be divided into two groups : the first 

indicating the motion of the wing in the tunnel and its reaction to the 
flowing gas and the second indicating the supplied power necessary to 
maintain the oscillation. A 1 4— channel recording oscillograph was used 
to record and correlate data from both groups with respect to time. 

The time history of the pitching motion of the wing was obtained 
from the signal of an electric strain gage which measured the torsio nal 
stress in- the torque rod. A light beam, reflected from a mirror 
fastened to the end of the wing and focused on a translucent scale, also 
gave an indication of the torsional amplitude of oscillation (fig. 1 ). 
Lift was obtained from the signal, of electric strain gages mounted on 
the bearing supports. 

The power required to maintain a fixed amplitude of oscillation 
was measured visually with a wattmeter which was incorporated in the 
electromagnetic— shaker circuit. At a correspon ding time, signals from 
an ammeter and voltmeter in the same circuit were recorded on the 
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oscillogram. At the point of wing resonance in pitching, the product 
of the volts and the amperes was equal to the measured watts. 


METHOD 


The oscillating lifts were measured by strain-sensitive elements 
on the hearing supports. The incidence of the wing was indicated by 
the shearing strain in the torsion spring. The phase relation of the 
lift to the incidence of the wing was obtained from the simultaneous 
recording of the lift and angular position on an oscillograph. The 
components of lift and the phase relation to wing incidence are 
illustrated in the following sketch. The frequency of oscillation is 
also obtained from the oscillogram. 



The aerodynamic damping moment, that is, the moment in phase with, 
the velocity, was obtained by the following procedures. S inc e a moment 
does work only with the component that is in- phase with the angular 
velocity, the energy required to sustain an oscillation is a measure 
of the dancing moment and is given by jtQ 0 a 0 sin cp 2 * Therefore the 

input energy, obtained from the wattmeter and corrected for tare nnd 
electric circuit losses, yielded the damping moment. To confirm this 
value, the aerody nami c damp ing moment and tare moment were obtained at 
maximum response (sin cpp = 1) from the ammeter reading and the physical 

and electric constants of the machine. This aerodynamic mome nt could 
then be corrected for the tare dampin g moment . Two other methods of 
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measuring damping in common use were considered. If either the 
logarithmic decrement or the damping from the response curve is used, 
the aerodynamic damping moment can he obtained. The results of these 
methods were in agreement. 

The aerodynamic moment in phase with the wing position may he 
obtained experimentally by observing the c hang es in the natural 
torsional frequency with Mach number. Such a method was used in 
reference 7 for an airfoil oscillating about its midchord. Since the 
model will oscillate in its "damped” natural frequency when the forcing 
stimulus is removed, the inphase moment may be determined (appendix A). 
With this knowledge, the shift of the inphase c ent er— of-pre s s ure 
position from the quarter chord may be determined by dividing by the 
inphase lift. 


ACCURACY OF EXPERIMENTAL RESULTS 


Before measuring the experimental oscillating lifts, the slope of 
the experimental static— lift curve of this wing was compared with the 
lift— curve slope for linear theory as demonstrated in figure 2. The 
theoretical dCjJda curve is the incompressible slope modified by the 


Glauert correction 



This corrected lift curve is comparable 


to K' sd for the case of oip = 0 as given by Dietze. Since there was 


some deviation of the experimental values from the theoretical values for 
the static case, the values of the experimental data for the oscillating 
case might also be expected to be less than the values indicated by the 
linear oscillatory theory. 


Quantitative measurements on high— frequency oscillations lead to 
certain problems not ordinarily encountered in static and low-frequency 
measurements. Deformation of the model under dynamic loading necessitates 
the correction of static calibrations. In the ca’se of torsional 
oscillations the static calibration varied from the dynamic conditions 
as the result of the difference between the static and dynamic deflection 
curves. This was corrected by applying a Holzer type of analysis to the 
dynamic configuration and relating the static strain relative to the 
midspan angle of incidence to the dynamic strain relative to the midspan 
angle of incidence. The deformation of the wing obtained from the Holzer 
analysis was approximately linear, the angle of attack of the side of the 
wing nearest the shakers being 9 6.5 percent of the incidence at the 
opposite end for wing— rod 1. The angle of incidence used in the 
calculation of the aerodynamic coefficients Is that at the midspan. 

In the measurement of the oscillatory lifts, the oscillograph 
records, particularly at the higher Mach numbers, were somewhat 
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distorted by wind-tunnel vibration and fluctuations in the flow of the 
air stream. To minimize the error due to this distortion, the oscillo- 
grams for both lift and angular position were read by taking an average 
over five consecutive cycles of an eight— point Fourier analysis. This 
analysis gave the amplitude of the lift vector and the magnitude of the 
angle of incidence as well as the phase between the lift and angular 
deflection. The precision of measurement as determined from repeated 
analysis from the same record gave lift values within 1 pound, the 
angle of incidence within 0.02°, and the phase angle within ±7°. Also, 
it was necessary to correct the phase relationship for any phase 
difference inherent in the galvanometer string. Unfortunately, the 
quadrature lift coefficient is very sensitive to variations of phase 
angle and therefore was subject to error. 

The correction due to dynamic magnification on wing-rod 1 was 
approximately 9 percent of the measured force (see appendix B). No 
attempt was made to correct the exper l m ental lifts and moments for the 
aerodynamic effect of this translation, since experiment and analysis 
showed the maximum displacement at midspan to be in the order of 
0.01 inch. This displacement for wing— rod 1 represents less than 
1 pound lift (approximately 1 percent total lift) under extreme tunnel 
conditions of high density and high Mach number. 

There was a possibility of the damping moment being somewhat in 
error due to the variations of the tare damping. Since the magnitude 
of the tare was from 25 to 60 percent of the measured value, small 
variations in the tare led to large variations in the aerodynamic— 
moment coefficients. 

A small random error was introduced by the necessity of making 
interpolations and extrapolations of the theory in order to compare the 
theory with the experimental results. The interpolation was made by 
cross— plotting the coefficients against Mach number for a given value 
of reduced frequency. The extrapolations were made by extending the 
curves of coefficients against reduced frequency given by Dietze and 
then cross— plotting the coefficients against Mach number at a given 
reduced frequency. This procedure tends to admit errors to the 
theoretical values. 


IMPROVEMENT OF TECHNIQUE 


For the future determination of oscillating air forces, there are 
various improvements to be effected. For more accurate results, 
dynamic calibrations for each specific test condition should be made. 
The difficulty in reading the oscillograms may be ameliorated by 
increasing the time scale or paper speed of the oscillograms. When 
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phase angles are “being measured, it has “been found that using the same 
type of signal, in similarly damped galvanometer elements will reduce 
any phase difference due to the galvanometers. It is also important to 
check that the response curves of these elements are n flat" for the 
frequency range being investigated. 

The “bending deflection of the wing— rod combination resulted in 
considerable dynamic magnification of the bearing reactions (see 
appendix B). The size of the correction is a function of the ratio of 
the forcing frequency to the natural bending frequency. To keep this 
magnification as low as possible, wings having very high bending 
frequencies compared to the forcing frequencies should be used. The 
tare damping in the Bystem should be a minimum, the evaluation being 
made at all test amplitudes . 


DISCUSSION AND RESULTS 


The theoretical coefficients of lifts and moments were calculated 
by using the theory of Dietze (references 3 ansi 4) . However, since the 
theoretical curves were determined to M = O .7 only, all reference to 
theory above that Mach number refers to an extrapolation of these 
curves. For comparison of the lift coefficients it is well to note 
that the inphase and quadrature components are functions of both the 
magnitude and the phase of the lift vector. “Under the test conditions 
it was impossible to treat these factors as independent variables. From 
figure 3 , it may be seen that a comparison can be made between the 
vector magnitudes of the experimental and theoretical lift coefficients. 
The difference between the experimental and theoretical values is shown 
by the short vector marked "difference.” It is shown that the 
experimental K’q-q is insensitive to variations from the theoretical 

phase angle cpp. Conversely, the K"gp is highly sensitive to such 

variations. Since the magn i tude of the lift but not the phase angle 
could be dete rmin ed accurately, the experimental, data which would be in 
closest agreement with theory would be magnitude of the total— lift 
vector. 

To show the range of values and agreement, the total theoretical 
lift has been plotted against experimental oscillating lift for wing- 
rod 1 (fig. *0 . The ratio of experimental total lift to theoretical 
lift was plotted against Mach number (fig. 5) • It is clearly shown 
that the theoietical values are greater than the experimental values. 

It may be noted that the order of deviation from the theoretical values 
is of the same magnitude as the deviation of the static— lift data shown 
in figure 2. In the same graph is plotted o^. showing the range of 

reduced frequencies in this set of tests. These and other data 
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pertaining to -wing-rod 1, that is, Mach number, coefficients, and so 
forth, are also shown in table I. 

Since the theory assuming an infinite wing in an incompressible 
flow with some form of compressibility correction is in common usage, 
it is of interest to compare the experimental, data with this theory. 

The ratio of experimental oscillating lift to this theoretical lift 
plotted against Mach number is given in figure 6. 

A comparison of the experimental, phase angle with theory given in 
figure 7 shows that the trends indicated by theory are followed. When 
it is remembered that sin is used in determining K"gj) (see fig. 3 ), 

it is only to be expected that the sign and order of magnitude of 
should deviate from the theory. Experimentally determined being 

insensitive to variations in the phase angle permits comparison with 
theory (fig. 8). It may be noted that the agreement is similar to that 
obtained for the magnitude of the lift vector. 

The opportunity was taken in this series of tests to examine the 
effects of small amplitudes of oscillation on the force coefficients. 

It has been suggested that for very small amplitudes some nonlinear 
effects associated with boundary layer may become significant. For the 
experiments herein reported no such effect was determined. The effect of 
various amplitudes of pitching oscillations is shown in figure 9» It can 
be seen that the ratio of coefficients for the four Mach numbers shown is 
relatively constant and independent of amplitude, indicating that the 
amplitude of oscillation had little effect on the lift coefficient. 

In view of the large corrections in the lift, angle of incidence, 
and phase angle necessary in reducing the data for wing— rod 3j it was 
deemed advisable to omit the presentation of data and Comparisons with 
theory. However, all trends were similar to those of wing-rod 1. 
Experiments were also conducted on wing— rod 2 which had a natural 
frequency of 45 cycles per second (midway between wing— rods 1 and 3)- 
The data obtained for this combination (not presented) are subject to 
question. It was found that the measured value of damping in pitch 
decreased with Mach number, approaching zero at M = 0.7. The ratio 
of the experimental to the theoretical oscillation lift coefficients 
also approached a minimum at M = .0.7. The data obtained for wing- 
rods 1 and 3 imply that this effect may be restricted to the particular 
test conditions for wing— rod 2. Unfortunately, time was not available 
to study this problem further. It might be noted that the phenomenon 
exhibits certain characteristics of a self— excited vibration. The 
reason for the behavior Is unexplained and It is hoped that future 
studies may answer the problem. 
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As previously mentioned , the experimental damping-moment coef- 
ficient in quadrature -with the torsional displacement was determined. 

A comparison of these experimental results with theory plotted against 
Mach number (fig. 10 ) shows the experimental values to he consistently 
low. This result has also "beep, noted hy Bratt and Chinneck (refer- 
ence 7 ) for subsonic speeds. Some of the variation in the experl mental 
coefficients may he attributed to fluctuations of the tare damping. 

The inphase moment coefficients E'pp were determined as shown in 

appendix A. Since, in the present case, the airfoil was oscillated 
about its quarter chord, the inphase moment was small. Thus, the 
change in natural torsional frequency may have been within the accuracy 
of the experiment. Therefore, a comparison for K’^ of experiment 

with theory shows little agreement (table I). However, for wing— rod 1 
which had the weakest spring and consequently the greatest sensitivity 
to aerodynamic moments, at the high Mach numbers (0.7 to 0 . 8 ) a 
negative oscillating moment occurred on the airfoil. This would indicate 
a rearward shift of the inphase center— of-pressure position of about 
3 percent chord. This center-of-pressure shift was predicted 
theoretically, but to a lesser degree (fig. 11) . 


coNCLmma remarks 


It is believed that the method of obtaining oscillating lifts and 
moments due to pitch as described in this paper can be refined to give 
sufficient accuracy for the experimental determination of those 
quantities in the transonic region. 

Good agreement of the experimental magnitude of the oscillating 
lift vector was obtained with the magnitude of the theoretical lift 
vectors of Dietze. Similar agreement was obtained between the 
experimental coefficients in phase with the wing displacement and the 
corresponding theoretical values. In general, both experimental 
coefficients were lower than the corresponding theoretical values. 

The trend with Mach number and reduced frequency indicated by 
theory for the phase between lift and incidence was obtained 
experimentally. The sensitivity of the quadrature components of lift 
to the variation in the experimental phase angle was such that the 
values obtained were questionable. 

The amplitude of torsional oscillation apparently has little effect 
on the lift coefficients in' the range of O.OOlt- to 0.014- radian. 
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A comparison of the experimental damping-moment coefficients with 
the theoretical values showed the experimental values to he consistently 
low. While quantitative agreement of the inphase moment coefficients 
with the theoretical values is not obtained, the movement of the inphase 
center of pressure tends to follow the behavior indicated by theory. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , November 30, 19^9- 
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AFPENDH A 

DETEEMEKATION OF THE INPHASE MOMENT AND 
THE CENTER— OF— PRESSURE SHUT 


The dete rminat ion of the inphase moment , when the change of natural 
frequency of oscillation of a wing is known, is as follows: 

The definition of the moment is the air spring times the wing 
angular deflection, namely AK^a 0 . The air spring is determined as 

follows: The natural frequency is 


“L = 



the frequency with air spring force is 


<!>> = 



thus 


_ % 

< 02 / ' s 6 + 


then 


&0 



(Al) 


Now, since the air spring is known, the inphase moment is easily 
calculated. 

From inspection of equation (Al), when AEg is positive, a>> 
and. there is a rearward shift of the center of pressure. When AKq is 
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negative, cop > o>i and there is a forward shift of the center of 

pressure. The position of the center of pressure referred to the 
quarter chord is dete rmin ed by the ratio of the inphase moment to the 
inphase lift force, namely 


C.P. 


AKeop 
2L 0 cos cpq 


This gives the center— of— pressure shift from quarter-chord position in 
percent chord. - 

For comparison with theory, the center-of— pressure shift was 
calculated by the theory of Dietze (reference 3) as follows: 


C.P. 


Q 


2L q cos cpq 

it pY 1 2 0 „ '( 3cv C 

- - — AZh,l-=f- 

2 2 °\ 8 


- K* 


dd 


y 2 /o^ 2 

U l AZOo^— - 

1 I *** ~ k, dd 

2 1 m 2 _ rrl 
2®r K SD 


Kt ED ~ 0, 373 
^'SD " “V 2 
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appendix: b 

DYNAMIC MAGNIEECATION OF LIFT FORCE 


When an elastic tody is subjected to an oscillating force, it is 
possible to determine the amplitude response of the resulting vibration. 
In the case of a simply supported undamped beam carrying a uniformly 
distributed oscillating load (reference 8) it may be shown that 


where 


7 = 7 b 



i=l 


sin a in ^ 
L 

$ - &)! 


= Mol* 

it^EI 


(Bl) 


L*> 


and y a is approximately equal to the static midspan deflection under a » 

uniform load ¥ 0 l. In view of this relation, and recognizing that the 

deflection curve is a function of the loading, it follows that the 
reactions at the supports are also dependent upon the deflection curve. 

This fact may be shown as follows for a simply supported beam having an 
arbitrary loading: 


7 = 



dx dx 


d^y 

dx^ 


<i) 


_L M 

EI dx 



dx 


dx 
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d 3 v = J M = JL 

3 El dx EE 
dx° 


where V, the vertical shear, is given by E 
span loading. 



dz and A(x) 


is 


It may he shown that 


lim d^y _ _B_ 
x-^o ^3 - H 


(E2) 


where E is the reaction at the support. Using the deflection curve 
from equation (Bl) , the reactions for a simply supported, undamped heam 
subjected to a uniformly distributed oscillating load are given by 



Treatment of this nature roughly approximated the physical 
conditions of the experiment and gives an indication of the order of 
the correction. For the higher frequency combinations, a closer approach 
to actual conditions is obtained by considering the wing— rod as a two- 
span continuous beam built in at the wall. The reactions of this 
structure were determined by using the method of moment distribution. 

By assuming a given distributed load the deflection curve was obtained. 
The inertia effects m afij were then known for that loading. Using the 
combination of the distributed and inertia loads a new deflection curve 
was obtained. By use of successive approximations the total loading was 
obtained. The ratio of the applied (aerodynamic) load to the measured 
load is then known. 

For the low— frequency wing— rod combination, this analysis, the 
results from equation (B 3 ), and experimentally determined correction 
factors were in agreement. The variation in the results of these 
methods for wing— rod Z, having a frequency ratio od/cuj, of O.65, was 
such that only the order of magnitude of the correction factor could 
bo determined. 
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Figure 1 .— 


Diagrammatic view of wing in the Langley 4 ^-— foot flutter 


re search tunnel. 
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Figure 2.— Slope of the static— lift curve as a function of Mach number 

for wing— rod 1. 
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Figure 3 .— against shewing oan^arisons of experimental lift coefficients and phase angle 

with theory. P = 0.00393* ° 0.013 for vingHrod 1. 
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Experimental lift, lb 


Figure 4 .— Theoretical lift against experimental lift far wing-rod 1 . 
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Mach number 


Figure 5.- Ratio of experimental lift to theoretical lift against Mach 
number and reduced frequency for ving-rod 1. 
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Theoretioal incompressible lift 
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Figure 6.- Ratio of experimental lift to theoretical lift for 
incompressible flow at various Mach numbers and reduced 
frequencies for wing-rod 1. 
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Figure 7.— Theoretical and experimental phase-angle relationship to 
Maoh number for ving-rod 1 . “ 1?4 radians per second. 
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Pitching displacement, radians 

Figure 9*— Ratio of experimental K*gp to theoretical K*gj) against 
pitching displacement for wing-rod 1. 
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Figure 10.— Ratio of experimental to theoretical K ,, - DD against 

Mach number and reduced frequency far wing-rod 1. 









Haoh number 


Figure 11.— Exper ime ntal and theoretical center— of— pres sure shift from, 
the 25-percent— chard position against Maoh number and reduced 
frequency far wing-rod 1. 





